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When a metal undergoes a continuous quantum phase transition, non-Fermi 
liquid behaviour arises near the critical point. It is standard to assume that all 
low-energy degrees of freedom induced by quantum criticality are spatially ex- 
tended, corresponding to long-wavelength fluctuations of the order parameter. 
However, this picture has been contradicted by recent experiments on a proto- 
type system: heavy fermion metals at a zero-temperature magnetic transition. 
In particular, neutron scattering from CeCug xAu^; has revealed anomalous dy- 
namics at atomic length scales, leading to much debate as to the fate of the local 
moments in the quantum-critical regime. Here we report our theoretical finding 
of a locally critical quantum phase transition in a model of heavy fermions. The 
dynamics at the critical point are in agreement with experiment. We also argue 
that local criticality is a phenomenon of general relevance to strongly correlated 
metals, including doped Mott insulators. 

Quantum (zero-temperature) phase transitions are ubiquitous in strongly correlated met- 
als; for a recent review, see ref. |l|. The extensive current interest in metals close to a second- 
order quantum phase transition has stemmed largely from studies of high-temperature su- 
perconductors. In these systems, however, it has been hard to actually locate the putative 
quantum critical points (QCPs). The situation appears to be simpler in some related families 
of strongly correlated metals. In particular, there are many heavy fermion metals which can 
be tuned between an antiferromagnetic (AF) metal and a paramagnetic metal. In recent 
years QCPs have been explicitly identified in a number of stoichiometric or nearly stoi- 



chiometric systems. Non-Fermi liquid behavioursi are observed, usually through transport 
and thermodynamic measurements, in the quantum critical regime. ExampleslHllI include 
CePd2Si2, CeCue-xAuj,, CeNi2Ge2, and YbRh2Si2. 

The most fundamental question regarding any condensed matter system concerns the 
nature of their low-energy excitations. The traditional theory for metals close to a QCP 
assumes that the only important degrees of freedom are long-wavelength fluctuations of the 
order parameter. For a magnetic QCP, such long-wavelength fluctuations are usually called 
paramagnons. This theory has its origins back in the 1960s in the description of liquid 



helium-3; for a historical perspective, see ref. 12. It was formulated in the modern lan- 



guage of critical phenomena in refs. |T^ and |T4[ The assumption that the long-wavelength 
paramagnons are the only critical degrees of freedom, adopted from the theory of classi- 
cal (finite-temperature) phase transitions, turns out to have three important consequences. 
Firstly, the dynamical spin susceptibility at generic wavevectors (wavevectors far from the 
AF ordering wavevector) has the usual Fermi-liquid form, i.e., it is linear in frequency. Sec- 
ondly, even close to the ordering wavevector, the dynamical spin susceptibility is also linear 
in frequency. Thirdly, there is a violation of the so-called uj/T scaling (defined in the next 
section) . 

This picture has been subjected to experimental test in recent years. Heavy fermion 
metals have undergone by far the most systematic study. A major puzzle has been raised by 
neutron-scattering experimentsiHl^ on CeCug-ajAua;, which goes from a paramagnetic metal 
to an antiferromagnetic metal as x increases through a critical value, Xc ~ 0.1 (ref. ^). Close 
to the AF ordering wavevector, a fractional exponent (a < 1) appears in the frequency 
and temperature dependences of the dynamical spin susceptibility. The same exponent a 
describes the frequency and temperature dependences at generic wavevectors, over essentially 
the entire Brillouin zone. Finally, the dynamical spin susceptibility exhibits uj/T scaling. 
These experiments have lead to much debate as to the nature of the quantum critical points 
in heavy fermion metals, especially concerning the fate of the local magnetic moments in 
the quantum-critical regimeS^BIli^B. 
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Here we show that a dynamical spin susceptibihty of the form observed experimentally 
may arise at a locally critical quantum phase transition, where local critical degrees of free- 
dom co-exist with spatially extended ones. While we reach this conclusion by solving a 
microscopic model appropriate for heavy fermion metals, the possibility of such a novel 
QCP can be seen on general grounds. The universality class of a classical phase transition is 
entirely determined by statics^. The critical degrees of freedom must be spatially extended, 
as only long-wavelength fluctuations of the order parameter are energetically favourable. 
At a quantum phase transition, by contrast, both the statics and dynamics (i.e., quantum 
fluctuations) are important^. The effect of quantum fluctuations on local degrees of freedom 
is to introduce a coupling to certain dissipative baths. When this coupling is sufficiently 
strong, the local degrees of freedom can become critical. (A similar situation — local degrees 
of freedom driven critical by a dissipative bath — arises in macroscopic quantum tunneling 
problems^.) 



Locally critical quantum phase transition 

The microscopic model we consider is the Kondo lattice model0, illustrated in Fig. |l|a. 
At each lattice site, a local momentill interacts via an exchange coupling with the spin 
of any conduction electron sitting at the site. There are two important energy scales in 
the problemilii: The Kondo temperature Tk sets the scale below which an isolated local 
moment would be screened by the spins of the conduction electrons, while the RKKY in- 
teraction characterizes the induced coupling between two local moments. We will assume 
an average occupancy of less than one conduction electron per lattice site; as a result, all 
phases considered here are metallic. To systematically analyze the behaviour of the lo- 
cal degrees of freedom in this model, we apply the extended dynamical mean-field theory 
(EDMFT) developed in refs. p4|-p6|. As in the standard dynamical mean- field theory§, the 
correlation functions of the lattice problem in the EDMFT are calculated through a self- 
consistent impurity Kondo problem. The latter, illustrated in Fig. |l]b, describes a local 
moment simultaneously coupled to two dissipative baths. A fermionic bath accounts for all 
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temporal fluctuations arising from hopping of electrons between the local site and the rest of 
the lattice, while a bosonic bath represents the fluctuating magnetic fleld generated by the 
local moments at all other sites. The couplings to the two baths are Jk and respectively. 
Further details of the method are described later in the paper. 

We flnd two types of QCP, illustrated by the phase diagrams in Figs. ^ and ^ In each 
case, the tuning parameter 5 is the ratio of the RKKY interaction to the Kondo temperature. 
Increasing 5 has two effects. First, the magnetic correlations become more pronounced. The 
dynamical spin susceptibility x(Q,ct;) — where Q is the peak wavevector — diverges at some 
threshold value 5 = 6c, which deflnes a QCP separating a paramagnetic metal {6 < 6c) 
from an antiferromagnetic metal {6 > 6c). Second, the local Kondo physics governed by the 
effective impurity problem is also change Through self-consistency, an increase in 6 

causes an increase in the coupling g between the local moment and the fluctuating magnetic 
field. At some value 6 = 6^^^, the corresponding ratio g/Tx reaches a threshold value where 
the local Kondo problem becomes critical. The existence of such a critical point reflects two 
competing processes in the local problem (Fig. |T]b): quenching of the local moment through 
its Kondo coupling to the spins of the fermionic bath, and the tendency of the coupling g to 
polarize the local moment along the direction of the fluctuating magnetic fleld. At 5 = (5fo^, 
the fluctuating magnetic fleld just succeeds in preventing the spins of the fermionic bath 
from completely quenching the local moment, yielding a singular local susceptibility. This 
critical point is also marked by the vanishing of a local energy scale E^^^, reflecting a critical 
slowing down. E*^^ serves as an effective Fermi energy scale, in the sense that the local 
susceptibility has the usual Fermi liquid form (i.e., an imaginary part that is linear in u) 
only at energies below E*^^. 

At the flrst type of QCP, shown in Fig. 0, the lattice system orders before the effective 
local problem has a chance to become critical, i.e., 6c < 6^^^, so E*^^ remains flnite at the 
transition. The local moments are completely quenched at zero temperature over the entire 
paramagnetic region and also at the QCP. Their only vestiges are the Kondo resonances, 
which hybridize with the conduction electrons to form Landau quasiparticleJ^liiisli^. The 
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RKKY interaction manifests itself as an interaction between these quasiparticles, leading 
to a spin-density-wave (SDW) instability which drives the transitionll. As a result, the 
dynamical spin susceptibility has the same form as in the traditional theoryEi'0. We find 
this type of QCP when the spin fluctuations are three-dimensional. 

The second type of QCP is illustrated in Fig. ^ Here, 5c = Si^^, so the lattice system 
reaches its ordering transition precisely at the point where the local problem also becomes 
critical. Thus, E*^^ vanishes at the QCP, and two kinds of critical degrees of freedom co- 
exist: long- wavelength fluctuations of the order parameter, and local fluctuations originating 
from the local moments. The transition is locally critical. We find this type of QCP when 
spin fluctuations are two-dimensional. (In order for the ordering temperature T/v shown in 
Fig. 1^ to be nonzero for 6 > 6c, it is necessary to have an infinitesimal RKKY coupling in 
the third dimension.) 

At the locally critical point, the zero-temperature dynamical spin susceptibility has an 
anomalous frequency dependence at wavevectors q not only close to the ordering wavevector 
Q but everywhere else in the Brillouin zone as well: 

Here a is an anomalous exponent [an expression for which is given in equation ([T3| ) below]. 
For any wavevector q far away from Q, /(q) is non-zero and varies smoothly with q; for q 
close to Q, /(q) is proportional to (q — Q)^. 

At finite temperatures, 00°" in equation (|l]) is replaced by T"A^(u;/T), where A4{uj/T) 
is a scaling function whose form depends only on a. This has two experimentally testable 
consequences. First, the dynamical spin susceptibility at q = Q satisfies an u/T scaling: 

Second, the static uniform spin susceptibility has a modified Curie- Weiss form 

^(^) = Q + 

where G is a positive constant. 



Our results can be directly compared to neutron-scattering experiments on Au-doped 
CeCug-ajAuj^ at the critical concentration Xc ~ 0.1: (i) The experimental data have been 
fittedilS to the form of equation ([^), with an exponent a ~ 0.75. (ii) Within experimental 
resolution, /(q) goes to zero along lines in the three-dimensional Brillouin zone, implying 
that the magnetic fluctuations are two-dimensional in real spaceSSH. (iii) An uj/T scaling 
has been extensively reportedi0 (see also ref. ^). 

In addition, the modified Curie- Weiss form [equation (^] is knownil0 to fit the mag- 
netization data in CeCue-xAua;. This form also appears to describe the magnetization of 
some other heavy fermion metals near quantum criticalityil. One example@ is YbRh2Si2, 
which is already very close to being critical at ambient pressure and without any doping. 
The thermodynamic and transport properties! of this compound are very similar to those 
of CeCug-xAua; at X = Xc- 

From the local susceptibility we can also determine the NMR spin-lattice relaxation 
rate 1/Ti. Barring subtle form-factor cancellations, the relaxation rate is expected to be 
temperature-independent in the relevant temperature range: 

— ~ constant. (4) 

This prediction could be tested, for instance, through NMR measurements on the Cu sites 
in CeCug-ajAuj; or the Si sites in YbRh2Si2. 

Extended dynamical mean field theory 

We now outline the EDMFT analysis leading to the conclusions reported above. The 
Kondo lattice model illustrated in Fig. |l]a is specified by the Hamiltonian 

= tij C^ia Cjc + Jk Sj ■ S(, j + Jjj Sj ■ Sj, (5) 

ij,a i ij 

where s^^j is the conduction electron spin at site i. Within the equation (U) is 

mapped onto an effective single-site problem, illustrated in Fig. and represented by the 
Hamiltonian 

T-Cioc = Jk S ■ Sc + ^Ep Cp^ + g^S- ((f)p + (j)lp^ + X! ' ^p^ (^) 

p,a p p 



where Cpo- and describe the fermionic and bosonic dissipative baths, respectively. The 
self-consistent procedure goes as follows. We put in trial forms for Ep, Wp, and g, and 
solve the impurity Kondo problem to determine the electron self-energy S(co'), the "spin 
self-energy" M{lj), the local Green's function Gioduj), and the local susceptibility Xioc(^^)- 
(For definitions of these functions and further details of the calculational procedure, see the 
Methods section. We note here that the key approximation of the EDMFT is the assumption 
that the electron and spin self-energies are momentum- independent.) Self-consistency is 
imposed by demanding that Gioc and xioc are equal to the wavevector averages, respectively, 
of the lattice Green's function G(k, cu) and the lattice susceptibility 

^(q'"^ = MM + /, ' 

where /q is the Fourier transform of lij. For the purpose of determining the universal low- 
energy behaviour, we take the density of states of the fermionic bath near the Fermi energy 
to be a non-zero constant, and the spectral function of the fluctuating magnetic field to have 
a power-law dependence on frequency, with an exponent 7, for u below some cut-off scale A. 
The effective impurity Kondo problem is solved using a 1 — 7 expansionil'^lii. Two types 
of self-consistent solution are obtained for different forms of the "RKKY density of states" , 

P/(e)=E^(^-^q)- (8) 
q 

The first type of solution, illustrated by Fig. 0, arises when the RKKY density of states 
p/(e) increases from the lower band edge (at e = Jq) in a square-root fashion (i.e., as 



e — Jq ). At the QCP, the spin self-energy takes the form 

M{uj) = -iQ-iau, (9) 

where a is a positive, real constant. Through equation (|^), this yields a ^(q, ^) of the form 
given in the traditional theoryllllii. The linear frequency dependence of M{uj) reflects the 
physics that the spin fluctuations actually describe quasiparticle-quasihole pairs. 

The second type of solution, illustrated by Fig. ^, arises when pi{e) increases from zero 
at the lower band edge with a jump. The self-consistent local dynamical spin susceptibility 
at the QCP is singular, and is given by 



XiocM = ^ In^, (10) 



where 

1 

A = — — exp 

7rpo(/i) 



(11) 



Po(/^) being the bare conduction electron density of states at the chemical potential /i. The 
corresponding spin self-energy has the form 

M{iu) ^-Iq + A lu". (12) 

Here, A = (— 'i)°'AoA~", where Aq is defined by the condition that, for e G (/q,/q + Aq), 
p/(e) is approximately equal to its value at the lower edge, pi{Iq). The exponent a depends 
on both A and P/(/q): 

The spin self-energy acquires this anomalous frequency dependence because spin fluctuations 
can not only decay into particle-hole pairs, but can also couple to critical local modes. 
The dynamical spin susceptibility at zero temperature is then given by equation (|T|), with 
/(q) = Jq — Jq. At finite temperatures, the self-consistent solution gives a spin self-energy 

M{uj, T) ^ -/q + A T^'M (io/T) . (14) 

The scaling function is M.{uj/T) = (i27r)" exp[a 'ip{l/2 — ic<j/27rT)], where ip is the digamma 
function. This directly leads to equations (||) and (|), with 6 = /q=o — /q and B ^ 
AoA-"(27r)"exp[aV^(l/2)]. 

A square-root onset and a jump at the lower edge of the RKKY density of states are 
characteristic of magnetic fluctuations in three and two dimensions, respectively. These 
behaviours are exemplified by the RKKY density of states associated with nearest-neighbor 
coupling in a (three-dimensional) cubic lattice and its counterpart in a (two-dimensional) 
square lattice. While the details depend on the lattice type and the range of interactions, 
the square-root and jump onsets in the two cases are robust provided that Jq approaches Jq 
in a (q — Q)^ fashion. 
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Beyond the microscopic theory 

So far we have presented the resuhs of a specific microscopic approach: the extended 
dynamical mean field theory for the Kondo lattice model. We now turn to more general 
considerations concerning the existence and properties of locally critical quantum phase 
transitions. There are four main points that we wish to discuss. 

Firstly, even in the traditional theory there are hints, albeit very subtle, that two- 
dimensionality facilitates the realization of a locally critical point. In this theorylilQ, the 
zero-temperature spin fluctuations are given by the harmonic, or Gaussian, fluctuations of 
the order parameter. As a result, the zero-temperature, zero-frequency spin susceptibility 
is proportional to (q — Q)^^. In two dimensions, the corresponding local — i.e., wavevector- 
averaged — susceptibility is in fact logarithmically singular. (For two-dimensional fermions 
coupled to two-dimensional commensurate spin fluctuations, there are other kinds of singu- 
larities in the SDW description; see ref. ^.) This singularity can cause important non-linear 
effects in the dynamics of any non-trivial local degrees of freedom — such as local moments 
in heavy fermion metals — and can ultimately drive them critical. Our EDMFT analysis 
amounts to a microscopic prescription for these non-linear effects. We note that the situa- 
tion is very different for quantum critical points in insulators0'0. There the effective fleld 
theory is usually below its upper critical dimension, implying that the spin susceptibility is 
(q — Q)"^"^^, where the anomalous exponent rj takes a positive value. The corresponding 
local susceptibility is no longer singular in two dimensions. 

Secondly, the contrast between a locally critical transition and an SDW transition can 
also be seen at the level of a Ginzburg-Landau-like description. For an SDW transition, there 
is only one type of critical degrees of freedom, namely the long-wavelength fluctuations of 
the order parameter. The Ginzburg-Landau (GL) action is then written purely in terms of 
the frequency-dependent magnetization at wavevectors close to the ordering wavevector: 



This generalizes the standard 0^ theory for a classical phase transition in that the fluctuations 
take place not only in the D spatial dimensions, but also in time. The temporal fluctuations 




(15) 
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can be thought of as adding z dimensions, where z is the dynamical exponent: iSiw is a 0^ 
theory with an effective dimensionahty Dcfr = -D + 2; see refs. |T3|,|T^|I|. 

For a locally critical point, on the other hand, the GL action is no longer just a 
theory. There are extra critical modes which characterize the continuous disappearance of 
the Kondo resonance as the critical point is approached from the paramagnetic side. These 
modes are independent of the fluctuations of the order parameter. They are spatially local 
since the destruction of the Kondo resonance is a local phenomenon. The GL action will 
now contain, in addition to the 0^ component, 5iw, parts describing these new local critical 
degrees of freedom and their non-linear coupling to the long- wavelength modes. The explicit 
construction of these additional parts of the GL action is left for future work. 

Thirdly, consideration of a GL description makes it likely that the locally critical point 
is stable beyond our EDMFT approximation, provided that a < 1 (as is the case in 
CeCug-xAui: and YbRh2Si2). The crucial question is what happens when we allow an 
explicit q dependence in the spin self-energy. The dynamical exponent entering will be 
z = 2 /a. When a < 1, the effective dimension D^s = D + z = 2 + 2/a will be above 
the upper critical dimension of 4, so that all the non-linear couplings within will be 
irrelevant in the renormalizat ion-group sense0. As a result, the q-dependence of the spin 
self-energy from iSiw will be at most (q — Q)^. We also expect that the spin self-energy 
due to the coupling of the order-parameter field with the local modes will have a smooth 
q-dependence. As a result, the zero-temperature spin susceptibility should be proportional 
to (q — Q)^^ (i.e., the anomalous exponent 77 = 0). This implies that the singular local 
susceptibility and the concomitant destruction of the Kondo resonance at the critical point 
are robust in two dimensions. It should be noted that a rigorous proof of this stability awaits 
the explicit construction of the GL action for the locally critical point. 

Fourthly, GL considerations also provide a way to understand the 00 /T scaling [equa- 
tion (0)] at the locally critical point. At an SDW transition, the field theory, equation (^5|), 
is a 0^ theory with an effective dimensionality D^s = d+z. Whereas at QCPs in insulators00 
the dynamical exponent is usually 1, in a metallic environment z = 2, raising Dcs above 
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the upper critical dimension of 4 for the 0^ theory. All the non-linear couplings are then 
irrelevant, and uj/T scaling is absentSS. (A violation of uj/T scaling has been reportedEj 
near a QCP in Cei_xLaxRu2Si2.) 

At a locally critical point, the non-linear couplings among the long- wavelength modes 
remain irrelevant. However, the GL action also contains non-linear couplings that involve 
the additional local critical modes. These couplings are relevant, making the field theory 
an interacting one, therebyil allowing the uj/T scaling given in equation @. In terms of a 
suitably defined spin relaxation rate, the contribution from the 0^ component depends on 
temperature with a power greater than 1. On the other hand, the contribution from the 
coupling to the local modes is linear in temperature due to the "relevant" non-linear cou- 
plings. The total relaxation rate is then linear in temperature, as required^ ioi uj/T scaling. 
Again, the EDMFT analysis of the previous section amounts to a microscopic prescription 
for determining the relaxation rate and, more generally, the entire scaling function. 

Broader implications 

In addition to explaining the salient features of experiments on heavy fermion metals close 
to a quantum critical point, our results may have broader implications for other strongly 
correlated metals. The emergence of critical local modes in the case we have studied depends 
crucially on the formation of local moments. A local moment is nothing other than a 
strongly interacting /-electron orbital which cannot be doubly occupied due to a strong on- 
site Coulomb interactionill. It is formed at energies intermediate between the bare Coulomb 
interaction scale and the asymptotic low-energy limit. The existence of non-trivial local 
physics at such intermediate energy scales is in fact ubiquitous in strongly correlated metals, 
including metals close to a Mott insulator. This suggests the possibility of a locally critical 
point occurring in a doped Mott insulator. 

Methods 

We adopt the EDMFT as a conserving resummation of diagrams for finite dimensional 
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systems [Fig. 7(b) of ref. |2^. The spatial dimensionality enters, among other ways, via 
the form of the RKKY density of states defined in equation For any finite-dimensional 
system, the support to p/(e) is bounded and a stable paramagnetic solution exists. 

The EDMFT provides a self-consistent procedure. In the Bose-Fermi Kondo Hamilto- 
nian, equation (^), we insert trial forms for the spectra of both the vector-bosonic bath 
and fermionic bath, which can be parameterized in terms of the Weiss fields Xo^ Gq^, 
respectively: 



rE:^^--xA^i E:j^^-G„i., (16) 

The trial forms for the QCP are taken to be 

ImG'o(cu + iO+) = -ttA^o, 
Imxg ^(cu + iO^) = C\u!\'^sgnLj for \uj\ < A, 



(17) 



where Nq, 7, C, and A are parameters to be determined. 

Next, we solve the Bose-Fermi Kondo problem using a (1 — 7) expansion. This step 
determines the local spin susceptibility and the local conduction electron Green's function, 

Xioc(r) = -{TrS^{T)S^{0))ni,^, 

(18) 

Gioc(r) = -(T.c.(r)ct (0))„,„^, 



where r is the imaginary time. This step also specifies^ the spin self-energy M{lj), and the 
conduction electron self-energy S(co'): 

(19) 



Gloc{uj)' 



(The spin self-energy is defined in terms of an effective spin cumulant that is J-irreduciblecJ.) 

The (1 — 7) expansion we use follows the work of Smith and SiH, Senguptail, and Vojta 
et a/.S. To linear order in 1 — 7, there is a renormalization-group fixed point located at 
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Jk = C = {1 — 7)/2. In analogy to the standard application of the e expansionHll, we assume 
that the local spin susceptibility Xioc(t) takes a power-law form at the critical point. By 
setting Jk and C to their fixed point values, we can determine the exponent by matching 
the leading ln(r) term in the expansion of the assumed power-law form for Xioc(t) with the 
leading ln(r) term of the perturbative expression for Xioc(t), calculated from the Bose- Fermi 
Kondo Hamiltonian, equation (|]), to linear order in 1 — 7. This gives, for < 7 < 1, 

Xioc(^) = ^r(7) sin ^^^^^{-zu;)-\ (20) 

where F is the gamma function, and for 7 = 0, equation (p!0|). 

In addition, to linear order in 1 — 7, the conduction electron self-energy vanishes: 



0((l-7f). (21) 



Finally, we demand self-consistency, which amounts to the requirement that a local cor- 
relation function is equal to the wavevector average of the corresponding lattice correlation 
function. In terms of the RKKY density of states, p/(e) defined in equation (H), and the usual 
conduction electron density of states, poi^), the self-consistency condition can be written as 

(22) 

Gioc^ = / de — -. 

J u — e — L^cj) 

When pi{e) is proportional to — Iq near its lower band edge (e Jq), xioc is finite 
at the QCP. The local moments are not critical. 

The locally critical point arises when pi{e) has a jump at the lower band edge. Xioc 
is now singular at the QCP. The self-consistent solution for the parameters introduced in 
equations (|T7D is 

No = Po{p), 7 = 0+, C = 7rA, (23) 

with the expression for A being given in the main text, as are those for xioc and M. We also 
note that the relationship between 7 and a is indirect: it is specified by the self-consistency 
equation (|22|) along with equations ([T7|)-(p!9D. 
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The local susceptibility, given in equation ([T0|), is universal. The exponent a for the 
spin self-energy depends on the product Apj{Iq) [see equation (p^j. It can be seen from 



equation (11) that A is of the order of the Kondo energy. Since the QCP is reached through 



a competition between the RKKY and Kondo interactions, we expect that this product is 
close to unity, resulting in an a that is not too far from |. The precise value of a, however, 
is interaction-dependent: it depends on which point of the phase boundary in the RKKY 
interaction-Kondo energy parameter space is crossed as the system is tuned through the 
quantum phase transition. 

At finite temperatures, the self-consistent procedure outlined above can also be carried 
through, leading to the spin self-energy quoted in the main text. 
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FIGURES 

FIG. 1. A theoretical model of heavy fermions. a, The Kondo lattice model. At each lattice 
site i, a spin-^ local moment Sj interacts with a local conduction electron orbital Cjo- through an 
antiferromagnetic Kondo coupling of strength Jx- The hopping amplitude between the conduction 
electron orbitals at sites i and j is tij. The Fourier transform of tij gives the band dispersion ek, 
corresponding to which is a density of states po{e). The Kondo temperature Tk is the characteristic 
scale for the screening of an isolated local moment by the spins of the conduction electrons. For 
sufficiently small Jk, Tk ~ [/9o(a*)]~^ exp[— l//Jo(/^)'/x], where ^ is the chemical potential. The 
RKKY interaction between two local moments at sites i and j is /jj, the Fourier transform of 
which is Iq. b, The effective impurity Kondo model to which the Kondo lattice model is mapped 
in the EDMFT§ih@. The local moment is coupled to two dissipative baths, one fermionic and 
the other bosonic. The bosonic bath describes a fluctuating magnetic field. The couplings to the 
two baths are Jk and g, respectively. The energy dispersions of the two baths are determined 
self-consistently. 

FIG. 2. Schematics of a conventional quantum phase transition in Kondo lattices. The tuning 
parameter 5 is the ratio of the RKKY interaction to the Kondo energy scale Tk- When 5 is small, 
the local moments are completely screened by the spins of the conduction electrons, and the system 
is a paramagnetic metal. Increasing 8 strengthens the coupling between each local moment and the 
fluctuating magnetic field produced by all other local moments, as measured by the ratio q/Tk, 
where g parametrizes the effective local problem defined in Fig. ||b and equation @. At 5^^^, the 
corresponding q/Tk would reach the critical value for a phase transition in the local dissipative 
problem (see main text); at this point, a local energy scale E*^^ would vanish. Increasing 6 also 
increases the spin susceptibility at the peak wavevector, x(Q)) which diverges at a quantum critical 
point at 8 = 5c- Beyond 6c, there exists a finite transition temperature T/v below which the system 
is an antiferromagnetic metal. In this conventional case, 5c < 5f^^, i.e., E*^^ is finite at the critical 
point. The quantum-critical (QC) regime, denoted by the wavy line, is described by the traditional 
theory. 
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FIG. 3. Schematics of a locally critical quantum phase transition in Kondo lattices. The 
notation is the same as in Fig. |2[ The crucial new feature here is that 5c = S^^^, i.e., E*^^ vanishes 
at the quantum critical point. Local and spatially-extended critical degrees of freedom co-exist in 
the quantum critical (QC) regime. 
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Figure 1 (Si et al.) 
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Figure 2 (Si et al.) 




Figure 3 (Si et al.) 



